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The singular limit e — > of the S-matrix associated with the equation ied^(t)/dt 



H(t)ip{t) is considered, where the analytic generator H(t) € M n (C) is such that its spectrum 
q-i is real and non-degenerate for all t € R. Sufficient conditions allowing to compute asymptotic 

formulas for the exponentially small off-diagonal elements of the S"-matrix as e — > are 
explicited and a wide class of generators for which these conditions are verified is defined. 
These generators are obtained by means of generators whose spectrum exhibits eigenvalue 
crossings which are perturbed in such a way that these crossings turn to avoided crossings. 
The exponentially small asymptotic formulas which are derived are shown to be valid up 
to exponentially small relative error, by means of a joint application of the complex WKB 
method together with superasymptotic renormalization. The application of these results to 
the study of quantum adiabatic transitions in the time dependent Schrodinger equation and 
of the semiclassical scattering properties of the multichannel stationary Schrodinger equation 
closes this paper. The results presented here are a generalization to n-level systems, n > 2, 
of results previously known for 2-level systems only. 
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1 Introduction 



Several problems of mathematical physics lead to the study of the scattering properties of 
linear ordinary differential equations in a singular limit 



ieif/(t) = H(t)ip(t), teR, e^O, 



where the prime denotes the derivative with respect to t, ip(t) G C n , H(t) G M n (C), for all 
t. A system described by such an equation will be called a n-level system. Let us mention 
for example the study of the adiabatic limit of the time dependent Schrodinger equation or 
the semiclassical limit of the one-dimensional multichannel stationary Schrodinger equation 
at energies above the potential barriers, to which we will come back below. When the 
generator H(t) is well behaved at +00 and —00, the scattering properties of the problem can 
be described by means of a matrix naturally associated with equation ( |1.1| ), the so-called 
S-matrix. This matrix relates the behavior of the solution ip{t) as t — > —00 to that of ip(t) 
as t — > +00. Assuming that the spectrum o~(t) of H(t) is real and non-degenerate, 



a(t) = {ei(i) < ea(t) < ■ • ■ < e n {t)} G R, 
the S'-matrix is essentially given by the identity matrix 



n.2 



S = diag (sn(e), s 22 (e), • • • s nn {e)) + 0(e°°), where Sjj(e) = 1 + 0{e) as e -» 0, (1.3) 

provided H(t) is C°°, see e.g. |Fi| , [|w| - Moreover, if H(t) is assumed to be analytic, 

it was proven in various situations that the off-diagonal elements Sjk of S are exponentially 



decreasing @, [fj, jJKPj , |JP4 ] 



(1.4) 



as e — > 0. See also [[JP1 ], Q, Q, Q for corresponding results in infinite dimensional 
spaces. Since the physical information is often contained in these off-diagonal elements, it is 
of interest to be able to give an asymptotic formula for Sjk, rather then a mere estimate. 

For 2-level systems (or systems reducible to this case, see [ JP2 ], [[|], [ MN| ), the situation 
is now reasonably well understood, at least under generic circumstances. Indeed, a rigorous 
study of the S'-matrix associated with (1.1) when n = 2, under the hypotheses loosely stated 



above, is provided in the recent paper [JP4]. The treatment presented unifies in particular 
earlier results obtained either for the time dependent adiabatic schrodinger equation, see e.g 
[JP2] and references therein, or for the study of the above barrier reflexion in the semiclassical 



limit, see e.g. [FF],Q. Further references are provided in [ JP4 1 . As an intermediate result, 
the asymptotic formula 



s jk = g jk e- r > k/£ (1 + 0(e)) , e -» 0, 



(1.5) 



for j ^ k G {1, 2}, with G C and Re Y ^ > is proven in [JP4]. As is well known, to get 
an asymptotic formula for Sjk, one has to consider (|l.l| ) in the complex plane, in particular 
in the vicinity of the degeneracy points of the analytic continuations of eigenvalues e\(z) and 
e2(z). Provided the level lines of the multivalued function 



Im 



ei(z') — e2(z')dz' = est, 



(1.6) 



called Stokes lines, naturally associated with ( |1.1| ) behave properly in the complex plane, the 
so-called complex WKB method allows to prove (|1.5[). But more important, it is also shown 
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in [JP4] how to improve ( |1,5[ ) to an asymptotic formula accurate up to exponentially small 
relative error 



Sjk 



with g 



9jk( £ ) e 

= 5jfc + C(e) and 1/ 



(l + 0(c 



-k/s 



0, 



:i.7) 



e) = r j 7 s + C'(e 2 ). This is achieved by using a complex WKB 



analysis jointly with the recently developed super asymptotic theory jBcfl , [Rf, [JP2|. Note 
that when given a generator, the principal difficulty in justifying formulas ( |1.5|) and ( |1.7| ) is 
the verification that the corresponding Stokes lines ( |1-6| ) display the proper behavior globally 
in the complex plane, which may or may not be the case | JKP| . However, this condition is 
always satisfied when the complex eigenvalue degeneracy is close to the real axis, as shown 
in 0| . See also [ MN | and [|R| for recent related results. 

For n-level systems, with n > 3, the situation is by no means as well understood. There 
are some results obtained with particular generators. In Q, flCHl ], [|CH2|1 and |BE|| , certain 
elements of the S-matrix are computed if H(t) = H*(t) depends linearly on t, H(t) = A+tB, 
for some particular matrices A and B. The choices of A and B are such that all components 
of the solution ip(t) can be deduced from the first one and an exact integral representation 
of this first component can be obtained. The integral representation is analyzed by standard 
asymptotic techniques and this leads to results which are valid for any e > 0, as for the 
classical Landau-Zener generator. The study of the three-level problem when H(t) = H*(t) £ 
M^(R) is tackled in the closing section of the very interesting paper fHPfl . A non rigorous 
and essentially local discussion of the behavior of the level lines Im Jq ej(z') — ej.(z')dz', 
j ^ k = 1,2,3, is provided and justifies in very favorable cases an asymptotic formula for 
some elements of the S-matrix. See also the review |So], where a non rigorous study of 
fll.lD is made close to a complex degeneracy point of a group of eigenvalues by means of an 
exact solution to a model equation. However, no asymptotic formula for Sj^, j ^ k, can 
be found in the literature for general n-level systems, n > 3. This is due to the fact that 
the direct generalization of the method used successfully for 2-level systems may lead to 
seemingly inextricable difficulties for n = 3 already. Indeed, with three eigenvalues, one has 
to consider three sets of level lines Im Jg ej(z') — ek(z')dz' to deal with ( |1.1| ) in the complex 
plane, and the conditions they have to fulfill in order that the limit e — > can be controlled 



may be incompatible for a given generator, see [Fl|, [F2] and [HP]. It should be mentioned 
however, that there are specific examples in which this difficult problem can be mastered. 



Such a result was recently obtained in the semiclassical study [Ba] of a particular problem 
of resonances for which similar considerations in the complex plane are required. 

The goal of this paper is to provide some general insight on the asymptotic computation 
of the S'-matrix associated with n-level systems, n > 3, based on a generalization of the 
techniques which proved to be successful for 2-level systems. The content of this paper is 
twofold. On the one hand we set up a general framework in which asymptotic formulas for the 
exponentially small off-diagonal coefficients can be proven. On the other hand we actually 
prove such formulas for a wide class of n-level systems. In the first part of the paper, we give 
our definition of the S'-matrix associated with equation (|1.1|) and explicit the symmetries it 



inherits from the symmetries of H(t), for t € R (proposition 2.1). Then we turn to the de- 
termination of the analyticity properties of the eigenvalues and eigenvectors of H(z), z € C, 
which are at the root of the asymptotic formulas we derive later (lemma [O]) . The next step 
is the formulation of sufficient conditions adapted to the scattering situation we consider, 



under which a complex WKB analysis allows to prove a formula like (1.5) (proposition 4.1). 
The conditions stated are similar but not identical to those given in |IKP| or fHPfl . As a 
final step, we show how to improve the asymptotic formula Ql.5| ) to (1.7) by means of the 
superasymptotic machinery (proposition 5^ and lemma 5.2). Then we turn to the second 
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part of the paper, where we show that a wide class of generators fits into our framework 
and satisfies our conditions. These generators are obtained by perturbation of generators 
whose eigenvalues display degeneracies on the real axis, in the spirit of We prove that 
for these generators, in absence of any symmetry of the generator H(t), one element per 
column at least in the S'-matrix can be asymptotically computed (theorem |6.1| ), This is the 
main technical section of the paper. The major advantage of this construction is that it is 
sufficient to look at the behavior of the eigenvalues on the real axis to check if the conditions 
are satisfied. The closing section contains an application of our general results to the study 
of quantum adiabatic transitions in the time dependent Schrodinger equation and of the 
semiclassical scattering properties of the multichannel stationary Schrodinger equation. In 
particular, explicit use of the symmetries of the S-matrix is made to increase the number of 
its elements for which an asymptotic formula holds. In the latter application, further specific 
symmetry properties of the S'-matrix are derived (lemma 7.1). 



Acknowledgments: It is a great pleasure to thank Charles-Edouard Pfister for many en- 
lightening and fruitful discussions which took place in Marseille and Lausanne. The hospital- 
ity of the Institut de Physique Theorique de l'EPFL where part of this work was accomplished 
is acknowledged. 

2 Definition and properties of the S'-matrix. 

We consider the evolution equation 

ieif/(t) = H(t)i/>(t), t G R, e -» 0, (2.1) 

where the prime denotes the derivative with respect to t, ip(t) G C n , H(t) G M n (C), for all 
t. We make some assumptions on the generator H{t). The first one is the usual analyticity 
condition in this context. 
HI There exists a strip 

S a = {z G C| |Im z\ < a}, a > 0, (2.2) 

such that H{z) is analytic for all z G S a . 

Since we are studying scattering properties, we need sufficient decay at infinity. 

H2 There exist two nondegenerate matrices H{+),H{— ) G M n (C) and a > such that 

lim |t| 1+a sup \\H(t + is) - H(±)\\ < oo. (2.3) 

^i 00 \s\<a 

We finally give a condition which has to do with the physics behind the problem. 
H3 For t£R, the spectrum of H(t), denoted by o~(t), is real and non- degenerate 

a{t) = {ei(t) < e 2 (t) <■■■< e n {t)} C R (2.4) 

and there exists g > such that 

inf \ ej (t) - e k (t)\ >g. (2.5) 



As a consequence of H3, for each t E R, there exists a complete set of projectors Pj(t) = 
Pf(t) G M n (C), j = 1, 2, • • • , n such that 

n 

E^) sI ( 2 - 6 ) 

n 

H(t) = Y / e j (t)P 3 (t) (2.7) 
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and there exists a basis of C n of eigenvectors of H(t). We determine these eigenvectors (fj(t), 
j = 1, 2, • • • , n uniquely (up to a constant) by requiring them to satisfy 



H{t)<pj(t) = e j (t)i Pj (t) 
Pj{t)^{t) = Q, i = 1,2, 



, n. 



(2.8) 
(2.9) 



Explicitly if Xj{t)i J ' = 1) 2, ■ ■ • , n form a complete set of differentiable eigenvectors of H(t), 
the eigenvectors 



^(t) = e-So^'^'^t), s.t. ^(0) = ^.(0) 



with 



_ (^(t)\P 3 (t)^(t)) . 
||^(t)||= ' J 



, n 



(2.10) 



(2.11) 



verify (2.9). That this choice leads to an analytic set of eigenvectors close to the real axis 
will be proven below. We expand the solution ip(t) along the basis just constructed, thus 
defining unknown coefficients Cj(t), j = 1, 2, • • • ,rt to be determined, 



™ ft 



(2.12) 



The phases e _ i ->o £j ( t ^ dt I s (see H3) are introduced for convenience. By inserting ( 2,12j ) in 
(2A) we get the following differential equation for the Cj(t)'s 



k=l 



iA jk (t)/e 



Cfc(t) 



where 



and 



A jk(t) = / (e^O - e k {t'))dt' 



a jk (t) 



(<p j (t)\ p J<t)<Pk(t)) 

ym 2 



(2.13) 



(2.14) 



(2.15) 



Here (-|-) denotes the usual scalar product in C n . Our choice ( |2.9[ ) implies ajj(t) = 0. It is 
also shown below that the a.,fc(i)'s are analytic functions in a neighborhood of the real axis 
and that hypothesis H2 imply that they satisfy the estimate 



lim sup |t| 1+a \ajk(t)\ < oo. 



t— >±oo 



3¥* 



(2.16) 



As a consequence of this last property and of the fact that the eigenvalues are real by 
assumption, the following limits exist 



lim Cj(t) = Cj(±oo). 



t±oo 



(2.17) 



We are now able to define the associated S'-matrix, S £ M n (C), by the identity 

ci(-oo) 

C2( — OO) 

\c n (-oo)/ \c n (+oo)/ 



s 



/ci(-oo)\ /ci(+oo)\ 

C2( + Oo) 



(2.18) 
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Such a relation makes sense because of the linearity of the equation ( [2.13 ). It is a well known 
result that under our general hypotheses, the S-matrix satisfies 

S = I + 0(e). (2.19) 

Note that the j th column of the S-matrix is given by the solution of ( |2.13| ) at t = oo subjected 
to the initial conditions Cfc(— oo) = 5jk, k = 1, 2, • • • , n. 

In general, the S-matrix defined above has no particular properties besides that of being 
invertible. However, when the generator H{t) satisfies some symmetry properties, the same 
is true for S. As such properties are important in applications, we show below that if H (t) is 
self-adjoint with respect to some indefinite scalar product, then S is unitary with respect to 
another indefinite scalar product. Let J £ M n (C) be an invertible self-adjoint matrix. We 
define an indefinite metric on C n by means of the indefinite scalar product 

(.,.)j = (-|J->. (2.20) 

It is easy to check that the adjoint of a matrix A with respect to the (■, ■) j scalar product 
is given by 

A* = J~ 1 A*J. (2.21) 



Proposition 2.1 Let H{t) satisfy HI, H2 and possess n distinct eigenvalues Vt E R. Fur- 
ther assume H(t) is self-adjoint with respect to the scalar product (•, 

H{t) = H#{t) = J- l H*(t)J, VteR, (2.22) 

and the eigenvectors <f j(0) of H(0) satisfy 

(<P3(P),<Pj(P))j = Pj, Pj € {-1,1}, Vj = l,---,n. (2.23) 

Then the eigenvalues of H(t) are realMt € R and the S -matrix is unitary with respect to the 
scalar product (•, where R = R* = R^ 1 is the real diagonal matrix R = diag (pi, p%, ■ ■ ■ , p n ), 

S* = RS*R = S' 1 . (2.24) 
Remark: 

The condition ((fj(0), <fj(0))j = ±1 can always be satisfied by suitable renormalization pro- 
vided (0),<^-(0)) j t^O. 

The main interest of this proposition is that when the S'-matrix possesses symmetries, 
some of its elements can be deduced from resulting identities, without resorting to their 
actual computations. 

A simple proof of the proposition making use of notions discussed in the next section 
can be found in appendix. The above proposition can actually be used for the two main 
applications we deal with in section [?]. Note that in specific cases, further symmetry property 
can be derived for the S'-matrix, see section [/]. 



3 Analyticity properties 

The generator H{z) is analytic in S a , hence the solution of the linear equation (|2.l| ) ip{z) 
is analytic in S a as well. However, the eigenvalues and eigenprojectors of H(z) may have 
singularities in S a . Let us recall some basic properties, the proofs of which can be found 
in jKJ]. The eigenvalues and eigenprojectors of a matrix analytic in a region of the complex 
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Figure 1: The paths 0, 5 and 770 in S a \Q. 



plane have analytic continuations in that region with possible singularities located at points 
zq called exceptional points. In a neighborhood free of exceptional points, the eigenvalues are 
given by branches of analytic functions and their multiplicities are constant. One eigenvalue 
can therefore be analytically continued until it coincides at zq with one or several other 
eigenvalues. The set of such points defines the set of exceptional points. The eigenvalues 
may possess branching points at an exceptional zq, where they are continuous, whereas 
the eigenprojectors are also multivalued but diverge as z — ► zq. Hence, by hypothesis H3, 
the n distinct eigenvalues ej(t) defined on the real axis are analytic on the real axis and 
possess multivalued analytic continuations in S a , with possible branching points at the set 
of degeneracies O, given by 

Q = {zq\ ej(zo) = efc(zo),for some k,j and some analytic continuation} . (3.1) 

By assumption H2, $7 is finite, by H3, f2 n R = and f2 = Q, due to Schwarz's principle. 
Similarly, the eigenprojectors Pj(t) defined on the real axis are analytic on the real axis and 
possess multivalued analytic continuations in S a , with possible singularities at f2. To see 
more precisely what happens to these multivalued functions when we turn around a point 
zq £ ft, we consider the construction described in figure [l|. Let / be a multivalued analytic 
function in S a \fl. We denote by f(z) the analytic continuation of /(0) along some path 
j3 € S a \VL from to z. Then we perform the analytic continuation of f{z) along a negatively 
oriented loop 5 based at z around a unique point zq € O, and denote by f(z) the function 
we get when we come back at the starting point (if 8 is positively oriented, the construction 
is similar). For later purposes, we define rjo as the negatively oriented loop homotopic to the 
loop based at the origin encircling zq obtained by following (3 from to z, 6 from z back to 
z and (3 in the reverse sense from z back to the origin. We'll keep this notation in the rest 
of this section. It follows from the foregoing that if we perform the analytic continuation of 
the set of eigenvalues {ej(z)}" =1 , along a negatively oriented loop around zq G fi, we get the 
set {ej(z)}™ =1 with 

ej(z) = e ao{j) {z), j = 1, ■■■,n, (3.2) 

where 

a : {l,2,---,n}^{l,2,---,n} (3.3) 

is a permutation which depends on 770 • Similarly, and with the same notations, we get for 
the analytic continuations of the projectors around zq 

Pj{z) = P aoU) (z), j = !,•••, n. (3.4) 
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Let us consider now the eigenvectors l -fj{t). We define W(t) as the solution of 

n 

W'{t) = Y,Pj(t)Pj(t)W(t) (3-5) 

= K(t)W(t), W(0) = I, 
where t G R. It is well known [|Kj] ; [jKrfl , that W(i) satisfies the intertwining identity 

W(t)Pj(0) =Pj(t)W(t), j = l,2,---,n , Vt G R, (3.6) 
so that, if {(^j(0)}" =1 denotes a set of eigenvectors of H(0), the vectors defined by 

<p j (t) = W(t)<p j (0) (3.7) 

are eigenvectors of H(t). Moreover, using the identity Q(t)Q' (t)Q(t) = which is true for any 
differentiable projector, it is easily checked that condition ( |2.9| ) is satisfied by these vectors. 
The generator K(t) is analytic on the real axis and can be analytically continued in S a \£l. 
Actually, K(z) is single valued in S a \Q. Indeed, let us consider the analytic continuation of 
K(z) around zq G fi. We get from (|3~4[ ) that 



= Kotfi*), ( 3 - 8 ) 

so that 



3=1 3=1 



n 



= Y. P k^)Pk{z) = K{z). (3.9) 
fc=i 

Consequently, W(t) can be analytically continued in S a \£l, where it is multivalued and 
satisfies both ( |3.5|) and ( |3.6|) with z G S a \Q in place of t G R. Moreover, the relation 
between the analytic continuation VF(z) from to some point z G S a \£l and the analytic 
continuation W(z) is given by a monodromy matrix W(rjo) such that 

W(z) = W(z)W(t 1o ), (3.10) 

where t]q is the negatively oriented loop based at the origin which encircles zq G only, (see 
figure |l]). Note also that the analytic continuation VK(z) is invertible in S a \£l and W^ 1 (z) 
satisfies 

W^'iz) = -W~ 1 (z)K(z), W _1 (0) = I. (3.11) 



As a consequence, the eigenvectors ( |3.7|) possess multivalued analytic extensions in S a \Cl. 
Consider the relation 

H(z)L Pj (z)=e j (z) i p j (z), (3.12) 

obtained by analytic continuation from to some point z G S a \£l. When analytically 
continued along a negatively oriented loop around zq G ft, it yields 

H(z){pj(z) = H(z){pj(z) = ej(z)yj{z) = e CTo(j) (z)^(z). (3.13) 

Thus <fj(z) is proportional to (p ao (j\(z) and we introduce the quantity 9j(r]o) G C by the 
definition 

{p j (z) = e- ie iM<p aoU) (z), i = l,2,-.-,n. (3.14) 
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This is equivalent to (see ( 3.1 0| ) ) 



W( Vo )!p j (0)=e- ie ^ ( p ao{j) (0). (3.15) 
Let us consider the couplings ( 2.15[) . Using the definition (|3.7|) , the invertibility of W(t) and 



the identity (|3.6| ), it's not difficult to see that we can rewrite 

a jk (t)- , ten, (3.16) 

which is analytic on the real axis and can be analytically continued in S a \Q, where it 
is multivalued. Thus, the same is true for the coefficients c~(t) which satisfy the linear 
differential equation ( 2.13j ) and their analytic continuations satisfy the same equation with 



z G S a \£l in place of t G R. We now come to the main identity of this section, regarding the 
coefficients Cj(z). Let us denote by Cj(z) the analytic continuation of Cj(0) from to some 
z G S a \Q. We perform the analytic continuation of Cj(z) along a negatively oriented loop 
around zq G SI and denote by Cj(z) the function we get when we come back at the starting 
point z. 

Lemma 3.1 For any j = 1, • • • , n, we have 

c i (z)e- i 4»^ (tt)du/e e-^('»)=c ffoCj . ) (z) (3.17) 
where tjq, 9j(j]o) and o" (j) are defined as above. 
Proof: 

It follows from hypothesis HI that ip(z) is analytic in S a so that 

n r z n pi — 

Y,Cj(z)e- l Jo e ^)dule^.^ z) = Y^c j (z)e- l Jo 'iM**/*-^ = 

3=1 3=1 

^^^^e^^^^/^/o^oW)^/^-^^)^^^). (3.18) 

3=1 

We conclude by the fact that {ifj{z)}^ =1 is a basis. □ 
Remark: 

It is straightforward to generalize the study of the analytic continuations around one singular 
point of the functions given above to the case where the analytic continuations are performed 
around several singular points, since £1 is finite. The loop r]Q can be rewritten as a finite 
succession of individual loops encircling one point of only, so that the permutation o"o is 
given by the composition of a finite number of individual permutations. Thus the factors 
e -iOj{vo) i n (3.14) should be replaced by a product of such factors, each associated with one 



individual loop and the same is true for the factors exp(— i J ej(z)dz/e) in lemma 3J. This 



process is performed in the proof of theorem 6.1 



4 Complex WKB analysis 

This section is devoted to basic estimates on the coefficients Cj(z) in certain domains ex- 
tending to infinity in both the positive and negative directions inside the strip S a . We first 
consider what happens in neighborhoods of ±oo. It follows from assumption HI by a direct 
application of the Cauchy formula that 

lim sup \t\ 1+a \\H'(t + is)\\ < oo. (4.1) 
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Figure 2: The path of integration for Ajj t (z) (the x's denote points of Q). 

Hence the same is true for the single valued matrix K{z) 

lim sup \t\ 1+a \\K(t + is)\\ < oo. (4.2) 

Let < T G R be such that 

min Re z > —T, and max Re z < +T. (4.3) 

zett zett 

All quantities encountered so far are analytic in S a n {z||Re z\ > T}, and we denote by a 
"~" any analytic continuation in that set. As noticed earlier 

W'{z) = K(z)W(z), z e S a n{z\\Re z\ > T} (4.4) 

so that it follows from (M.2I) that the limits 



lim Wit + is) = W(±oo) (4.5) 

t— >±oo 



exist uniformly in s €] — a,a[. Consequently, see (|3.16), 



lim |t| 1+a sup \a jk (t + is)\ < oo, Vj, k e {1, • • • , n\. (4.6) 

*^ ±oc \s\<a 

Finally, for \t\ > T, we can write 

Im Ajk(t + is) = Im (^J ej(z)dz — J ek(z)dz^j 

+ jf Re (e CTj(j) (t + is') - e a . {k) (t + is'))ds', (4.7) 

where this equation is obtained by deforming the path of integration from to z = t + is 

into a loop rj based at the origin, which may encircle points of £1, followed by the real axis 
from to Re z and a vertical path from Re z to z, see figure §. Hence we have 

sup Im Ajfc(z) < oo, (4-8) 

zeS^nlzWRc z\>T} 

which, together with ( |4.6j ) yields the existence of the limits 

^ lim Cj (t + is) = Cj (±oo) (4.9) 

uniformly in s g] — a,a[. We now define the domains in which useful estimates can be 
obtained. 
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Figure 3: The path j3 along which the analytic continuation of Aj k (t) in Dj is taken. 



Definition: Let j E {1, • • • ,n} be fixed. A dissipative domain for the index j, Dj C S a \U, 
is such that 



sup Re z = oo, inf Re z = — oo 



zeD, 



and is defined by the property that for any z E Dj and any k E {1, 
7 fe C Dj parameterized by u e] — oo,t\ which links — oo to z 

lim Re ~f k (u) = -oo, -y k (t) = z 



u— >— oo 



with 



sup sup 

z&Dj mG]- oo,t] 



< OO 



(4.10) 

,n}, there exists a path 
(4.11) 

(4.12) 



and satisfies the monotonicity condition 

Im Ajfc(7 fe (ii)) zs a non decreasing function of u E] — oo, i]. 



(4.13) 



Suc/i a pai/i is a dissipative path for {jfc}. -ffere Aj k (z) is the analytic continuation of 



Ajk(t) = f\ej(t r ) - e k {t'))dt', t E R, 
Jo 



(4.14) 



in Dj along a path f3 described in figure [| going from to —T E R along the real axis and 
then vertically up or down until it reaches Dj, where T > is chosen as in (U-V - 
Remark: 

The finiteness of f2 insures the existence of such a path (3. 

Let c k (z), k = 1, 2, ■ • • , n, z E Dj, be the analytic continuations of c k (t) along the same 
path (3 which are solutions of the analytic continuation of ( 2,13| ) in Dj along (3 



»A fc i(z)/e 



ci(z). 



1=1 



We take as initial conditions in Dj 

lim c k (z) = lim c k {t) = 5 jk , k = 1, • 

Re z^-oo t^-oo 

and we define 

x k (z) = c k (z)e iK ^/ e , zeDj, k = !,•••, n. 



,n. 



(4.15) 

(4.16) 
(4.17) 
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Lemma 4.1 In a dissipative domain for the index j we get the estimates 



sup \xj(z) - 1| = 0(e) (4.18) 

sup |x fc (z)| = 0{e), Vk^j. (4.19) 

Remark: 

The real axis is a dissipative domain for all indices. In this case we have Cj(t) = Cj{t). Hence 
we get from the application of the lemma for all indices successively that 

S = I + 0(e). (4.20) 

The estimates we are looking for are then just a direct corollary. 

Proposition 4.1 Assume there exists a dissipative domain Dj for the index j. Let rjj be a 
loop based at the origin which encircles all the degeneracies between the real axis and Dj and 
let o~j be the permutation of labels associated with rjj , in the spirit of the remark ending the 
previous section. The loop rjj is negatively, respectively positively, oriented if Dj is above, 



respectively below, the real axis. Then the solution of ( 2. IS ) subjected to the initial conditions 
Ck{— oo) = 5jt satisfies 

S . (j) (+oo) = e-^e- l ^ eAz)dz/£ (l + 0(e)) (4.21) 
c aj{k) (+oc) = O (ee lm 4 ^W^M SW (+»)~S^>(^)/ E j > (4 . 22) 

with hj £ [H7, Hj~], where Hj is the maximum, respectively minimum, imaginary part of 
the points at +oo in Dj 

H + = lim sup s, H~ = lim inf s. (4.23) 

*^+°° s \t+iseDj t^+oo s \t+iseDj 

Thus we see that it is possible to get the (exponentially small) asymptotic behavior of the 
element Sa(j),j of the S-matrix, provided there exists a dissipative domain for the index j. 
The difficult part of the problem is of course to prove the existence of such domains Dj, 
which do not necessarily exist, and to have enough of them to compute the asymptotic of 
the whole S'-matrix. This task is the equivalent for n-level systems to the study of the global 
behavior of the Stokes lines for 2-level systems. We postpone this aspect of the problem to 
the next section. Note that we also get from this result an exponential bound on the elements 
s tTj(k),j OI t ne S'-matrix, k ^ j, which may or may not be useful. If rjj encircles no point of 
Q, we cannot get the asymptotic behavior of s a .^j but we only get the exponential bounds. 
Since our main concern is asymptotic behaviors, we call the corresponding dissipative domain 
trivial. 
Remark: 

In contrast with the 2-level case, see [ [JP4 ], we have to work with dissipative domains instead 



of working with one dissipative path for all indices. Indeed, it is not difficult to convince 
oneself with specific 3-level cases that such a dissipative path may not exist, even when the 
eigenvalue degeneracies are close to the real axis. In return, we prove below the existence of 
dissipative domains in this situation. 
Proof: 

The asymptotic relation is a direct consequence of lemma , (|4.9[) , (|4.17| ) and the first part 
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of the lemma. The estimate is a consequence of the same equations, the second estimate of 
the lemma and the identity, for t > T, 



Im A.jk(t + is) 



Im 



ej(z)dz 



e k (z)dz 



'I, 



+ Jo Re ( e °*G)(* + is '^ ~ e ^( k )^ + is ')) ds> - 



(4.24) 



The path of integration from to z for Aj k (z) is deformed into the loop rjj followed by the 
real axis from to Re z and a vertical path from Re z to z. It remains to take the limit 
t -> +oo. □ 



Proof of lemma 4.1 



We rewrite equations ( 4.15| ) and ( 4.16| ) as an integral equation and perform an integration 
by parts on the exponentials 

aki(z) 



i=i 

n 

+ i£ £ 

Lm=l 



1=1 



ek(z - ei[z 



-e iAkliz)/£ ci(z) 



aki(z') 



e k (z') - ei(z') 



akl{z')ai m {z') jA km (z')/e 



oo e k (z') - ei(z') 



c m (z')dz'. 



(4.25) 



Since all eigenvalues are distinct in S a \£l, the denominators are always different from 0. Due 
to equation (|4.9| ), the height above or below the real axis at which we start the integration 
is irrelevant, so that we can use the symbol — oo as lower integration bound. Note that the 
integrated term vanishes at — oo. We have also used the identity 

A kl (z') + A lm (z') = A km (z'). (4.26) 

In terms of the functions x k we get, using the same identity, 



x k (z) 



aki(z) 



+ 

i=i 

n 

+ is 



l= l ek(z) - ei(z) 
a k i(z r ) 



xi(z) 



e k {z') - ei(z') 



e <*M*)-*M*')ye Xl ( z >)dz' 



i^J-oo e k (z') - ei(z') 
We introduce the quantity 

|||x|||j = SUp |x;(z)| 



a fe /(^)^m(4 ei (A, fe (,)-A, fc ( z 0)/ £xTO( / ) ^/ 



(4.27) 



(4.S 



and consider for each k the equation fl4,27j ) along the dissipative path j k (u) described in the 



definition of Dj, such that 



e i(A jfc ( 7 fc (t))-A jfc ( 7 fc (u)))/e 



< 1 



(4.29) 



when u < t along that path. Due to the integrability of the a k i(z) at infinity and the uniform 
boundedness of dj k (u)/du, we get the estimate 



\x k {z) - 5 kj \ < e\\\x\\\jA 



(4.30) 
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for some constant A uniform in z G Dj, hence 

HMHi < l + e|||a;|||jA (4.31) 
Consequently, for e small enough 

HWIIj<2. (4.32) 
And the result follows. □ 

5 Superasymptotic improvement 

All results above can be substantially improved by using the so-called superasymptotic renor- 
malization method [Be], Q, [ JP2| . The joint use of complex WKB analysis and superasymp- 



totic renormalization is very powerful, as demonstrated recently in [ [JP4| for 2-level systems, 
and allows, roughly speaking, to replace all remainders O(e) by 0(e _K//e ), where k > 0. We 
briefly show how to achieve this improvement in the case of n-level systems. 
Let H(z) satisfy HI, H2 and H3 in S a and let 

Sa = S a \ U r =l,... l p (J r U J r ), (5.1) 

where each J r is an open domain containing one point of ft, only in the open upper half 
plane. Hence, any analytic continuation Cj(z) of ej(t), t £ R, in S a is isolated in the 
spectrum of H(z) so that &j(z) is analytic and multivalued in S a , and the same is true for 
the corresponding analytic continuation Pj(z) of Pj(t), t £ R. Let a r be the permutation 
associated with the loop Q r based at the origin which encircles J r once, such that 

Zj(z) = e ar{j) (z), (5.2) 

with the convention of section |3[ The matrix K(z) is analytic and single valued in S a . 
Consider the single valued analytic matrix 

Hx(z,e) = H{z) - ieK(z), z e S a . (5.3) 

For e small enough, the spectrum of H\(z, e) is non degenerate Vz £ S a so that its eigenvalues 
ej(z,e) and eigenprojectors Pj{z,e) are multivalued analytic functions in S a . Moreover, for 
e small enough, the analytic continuations of e}j(z,e) and Pj(z,e) around J r satisfy 

e){z) = el rU) (z) (5.4) 
Pj{z) = Pl U) (z), (5.5) 
as can be easily deduced from (^^) by perturbation theory. Consequently the matrix 

m 

K 1 (z,£)=J2P}\z,e)P}(z,s) (5.6) 
i=i 

is analytic and single valued in S a - Defining the single valued matrix 

H 2 (z,e) =H(z) -ieK x {z,e), z € S a , (5.7) 

we can repeat the argument, for e small enough. By induction we set for any q € N, 

H q (z,e) = H(z)-i£K q ^(z,e) (5.8) 

m 

e) = ^if-^Ejif- 1 (*,£), ZG§ a (5.9) 
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for e is small enough. We have 

rn 

H q (z,e) = ^e%z,e)Pf(z,e), (5.10) 

3=1 

where the eigenvalues and eigenprojections are multivalued in S a and satisfy 

e^(z,e) = el (j) (z,e) (5.11) 

P q (z,e) = Pl u) (z,e), j = l,---,n, (5.12) 

with the notations of ( |5.2[ ). We quote from [JP4], [JP2] the main proposition regarding this 
construction. 

Proposition 5.1 Let H(z) satisfy HI, H2 and H3 in S a and let S a be defined as above. Then 
there exist constants c > 0, e* > and a real function b(t) with lim^-i-oo \t\ 1+a b(t) < oo, 
such that 

\\K q (z,e) - K q ^{z,e)\\ < b(Re z)e q c q q\ (5.13) 
\\K q (z,e)\\ <b(Re z), (5.14) 

for all z € S a , all e < e* and all q < q*(e) = [1/ece], where [y] denotes the integer part of y 
and e is the basis of the neperian logarithm. 

We can deduce from this that in S a 

e)(z,e) = ej(z) + 0{e 2 b{Re z)) (5.15) 
P]{z,e) = e j (z) + 0(e(Re z)), \/q<q*(e). (5.16) 

We introduce the notation f q *^ = f* for any quantity f q depending on the index q and 
we drop from now on the e in the arguments of the functions we encounter. We define the 
multivalued analytic matrix W*(z) for z G S a by 

Wj{z) = K*(z)W*(z), W»(0) = I. (5.17) 



Due to the above observations and proposition |5.l| , W*(z) enjoys all properties W(z) does, 
such as 

W*(z)P*(z) = P*(0)W*(z) (5.18) 

W*(z) = W*(z)W*(C r ) (5.19) 
and, uniformly in s, 

lim WJt + is) = WJ oo). (5.20) 

Thus we define for any z G S a as set of eigenvectors of H*(z) by 

<p* j (z)=W.(z)<p* j (0), (5-21) 

where 

iT*(0)^*(0) = e*(0)^(0), j = l,---,n, (5.22) 
and which satisfy 

^(0)=e-^^V;. a) (0), (5-23) 
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with 



e](( r ) = 9(( r ) + o(e) e c. 



Let us expand the solution of (2.1) on this multivalued set of eigenvectors as 



n 



(5.24) 



(5.25) 



Since the analyticity properties of the eigenvectors and eigenvalues of H*(z) are the same as 



those enjoyed by the eigenvectors and eigenvalues of H(z), we get as in lemma 3.1 

-*/ \ -* f, e*Au)du/e -iO*(r r ) * , n w ^- c 

Cj(z)e J e 3 ^= c ^ (i) (^), \/zeS a . 



(5.26) 



Substituting ( 5.25| ) in (|2.1|), we see that the multivalued coefficients c*(z) satisfy in S a the 
differential equation 



k=i 



where 



A* jk (z) = I e*(z')-et(z')dz> 



and 



(^(^(O)!^^)^)^^)- 1 ^,!^)-^^))^^)^^)) 

ll^(0)|| 2 



(5.27) 



(5.28) 



(5.29) 



to be compared with ( 3.16| ). The key point of this construction is that it follows from 

(5.30) 



proposition 5A with q = q*(e) that 

\a* k {z)\ < 26(Re z)e~ K/£ , Vz G S a 
where k = 1/ec > 0, and from (|5.15|) that 

Im A* k (z) = Im A jk (z) + 0(e 2 ), 
uniformly in z G S a . Thus, we deduce from ( |5.30| ) that the limits 
lim c*(t + is) = c*(±oo), j = l,---,n, 

t^±oo J J 



(5.31) 



(5.32) 



exist for any analytic continuation in S a - Moreover, along any dissipative path j k (u) for 
{jk}, as defined above, we get from (|5.31|) 



e i(A* ft ( 7 fe(t))-A* A ( 7 fc(«)))/ £ 



0(1) 



(5.33) 



so that, reproducing the proof of lemma 4.1 we have 

Lemma 5.1 In a dissipative domain Dj, i/c^(— oo) = c k (— oo) = <5fcj, t/ien 
c f j {z) = l + 0(e- K ' E ) 
<&">W%(z) = O (e- K / £ ) , Vk ± j, 

uniformly in z £ S a . 



(5.34) 
(5.35) 
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This lemma yields the improved version of our main result. 



Proposition 5.2 Under the conditions of proposition [fj and with the same notations. If 
c* k {—co) = djk, then 



c; j(j) (+oo) = e-^e- l ^^ Z)dz/£ (l + 0(e-^)) (5.36) 

4 jW (+o°) = o( e -"/«e to 4^ W * / ^ ( ^« ( ^ ) ^« ( ^ ,)/e V (5.37) 

Note that we may or may not replace ej(z) by e*j(z) in the estimate without altering the 
result. It remains to make the link between the 5-matrix and the c£(+oo)'s of the proposition 
explicit. We define /3j by the relations (H*(z) and H{z) coincide at ±oo), 

tp* j (±oo) = e~ ip i ± (p j (±oo). (5.38) 
By comparison of flOD and (|2~rp we deduce the lemma 

Lemma 5.2 If c k {t) and cj£(i) satisfy c k {— oo) = c£(— oo) = <5jfc, i/ien, i/ie element kj of the 
S-matrix is given by 

S kj = C k (+00) = e -^r-/ 3 r) e - i /o + ^^( i Ve fc (0^7 £e -/L^(0- J ( t 'K'/ £4(+oo) 

= e-^c^+oo), (5.39) 
TOi/i 0*+ = 0(e) and ej(tf) - ej (t')dt'/e = 0(e), i.e. e" ia ^ = 1 + 0(e). 
Remarks: 



i) Proposition 5^ together with lemma [5^ are the main results of the first part of this paper. 

ii) As a direct consequence of these estimates on the real axis we have 



and 



s jk = (e-«/ £ ) , Vk^j, (5.40) 
Sjj = e ~ ia h (l + O (e- K ' £ )) . (5.41) 



iii) It should be clear from the analysis just performed that all results obtained hold if the 
generator H(z) in fl2.1|) is replaced by 



H(z, e) = H Q {z) + 0(e6(Re z)), (5.42) 
with b(t) = 0(l/t 1+a ), provided Hq(z) satisfies the hypotheses we assumed. 



6 Avoided crossings 

We now come to the second part of the paper in which we prove asymptotic formulas for 
the off-diagonal elements of the S-matrix, by means of the general set up presented above. 
To start with, we define a class of n-level systems for which we can prove the existence of 
one non trivial dissipative domains for all indices. They are obtained by means of systems 
exhibiting degeneracies of eigenvalues on the real axis, hereafter called real crossings, which 
we perturb in such a way that these degeneracies are lifted and turn to avoided crossings 
on the real axis. When the perturbation is small enough, this process moves the eigenvalue 
degeneracies off the real axis but they remain close to the place where the real crossings 
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occurred. This method was used successfully in 0| to deal with 2-level systems. We do not 
attempt to list all cases in which dissipative domains can be constructed by means of this 
technique but rather present a wide class of examples which are relevant in the theory of 
quantum adiabatic transitions and in the theory of multichannel semiclassical scattering, as 
described below. 

Let H(t,5) G M n (C) satisfy the following assumptions. 

H4 For each fixed S G [0, d], the matrix H(t,5) satisfies HI in a strip S a independent of 5 
and H(z,5), d/dzH(z,5) are continuous as a functions of two variables (z, 5) G S Q x [0,d]. 
Moreover, it satisfies H2 uniformly in 5 G [0,d], with limiting values H(±,5) which are 
continuous functions of 5 G [0,d]. 

H5 For each t G R and each 5 G [0,d], the spectrum of H(t,5), denoted by a(t,5), consists 
in n real eigenvalues 

a(t,S) = {e 1 (t,5),e 2 (t,S),---,e n (t,S)} CR (6.1) 

which are distinct when 5 > 

ei(t,S) < e 2 (t,5) <■ ■ ■< e n (t,5). (6.2) 

When 5 = 0, the functions ej(t,0) are analytic on the real axis and there exists a finite set 
of crossing points {t\ < t 2 < • • • < t p } G R, p > 0, such that 

i) Vt < h, 

ei(t,0) < e 2 (t,0) <■■■< e n {t,0). (6.3) 

ii) Vj < k G {1, 2, • • • , n}, there exists at most one t r with 

e i (t P ,0)-e fc (t rj 0) = > (6.4) 
and if such a t r exists we have 

^(e i (t P> 0)-e fe (t r ,0))>0. (6.5) 

Vj G {1, 2, • • • , n}, the eigenvalue ej(t, 0) crosses eigenvalues whose indices are all superior 
to j or all inferior to j. 
Remarks: 

i) The parameter 5 can be understood as a coupling constant controlling the strength of the 
perturbation. 

ii) The eigenvalues ej(t,0) are assumed to be analytic on the real axis, because of the degen- 
eracies on the real axis. However, if H (t, 5) is self adjoint for any 5 G [0, d], this follows from 
a theorem of Rellich, see Q. 

iii) We give in figure |] an example of pattern of crossings with the corresponding pattern of 
avoided crossings for which the above conditions are fulfilled. 

iv) The crossings are assumed to be generic in the sense that the derivative of e» — are 
non zero at the crossing t r . 

v) The crossing points {ti,t 2 , ■ • • ,t p } need not be distinct, which is important when the 
eigenvalues possess symmetries. However, for each j = 1, • • • , n, the eigenvalue ej(t,S) ex- 
periences avoided crossings with ej + i(t,5) and/or ej-i(t,5) at a subset of distinct points 
{t n , ■ ■ ■ , t rj } C {t\, t 2 , ■ ■■ ,t p }. 

We now state the main lemma of this section regarding the analyticity properties of the 
perturbed levels and the existence of dissipative domains for all indices in this perturbative 
context. 
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Figure 4: A pattern of eigenvalue crossings (bold curves) with the corresponding pattern of avoided 
crossings (fine curves) satisfying H5. 



Lemma 6.1 Let H(t,5) satisfy H4 and H5. We can choose a > small enough so that the 
following assertions are true for sufficiently small 5 > 0: 

i) Let {t ri , ■ ■ ■ ,t r } be the set of avoided crossing points experienced by &j(t, 8), j = 1, • • • , n. 
For each j, there exists a set of distinct domains J r G S a , where r G {n, • • • ,rj}, 

J r = { z = t + is\0 < \t - t r \ < L, < g < s < a'}, (6.6) 

with L small enough, a' < a and g > such that e, (— oo,<5) can be analytically continued in 

S J Q = S a \ U r=ri ,..., r . (J r U Jr). (6.7) 

ii) Let t r be an avoided crossing point of ej(t,S) with ek(t,5), k = j ± 1. Then the analytic 
continuation of ej(t r ,5) along a loop based at t r G R which encircles J r once yields ej(t r ,5) 
back at t r with 

ej(t r ,5) = e k (t r ,5). (6.8) 

Hi) For each j = 1, • • • , n, there exists a dissipative domain Dj above or below the real axis in 
S a Pi{z = t + is\ \s\ > a'}. The permutation aj associated with these dissipative domains (see 
proposition |^. j\ ) are all given by aj = a where a is the permutation which maps the index of 
the k til eigenvalue ej(oo, 0) numbered from the lowest one on k, for all k G {1, 2, • • • , n}. 

Remarks: 

i) In part ii) the same result is true along a loop encircling J r . 

ii) The dissipative domains Dj of part iii) are located above (respectively below) all the sets 
J r (resp. J r ), r = 1, • • • ,p. 

iii) The main interest of this lemma is that the sufficient conditions required for the exis- 
tence of dissipative domains in the complex plane can be deduced from the behavior of the 
eigenvalues on the real axis. 

iv) We emphasize that more general types of avoided crossings than those described in H5 
may lead to the existence of dissipative domains for certain indices but we want to get dis- 
sipative domains for all indices. For example, if part iii) of H5 is satisfied for certain indices 
only, then part iii) of lemma O is satisfied for those indices only. 
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v) Note also that there are patterns of eigenvalue crossings for which there exist no dissipa- 
tive domain for some indices. For example, if ej(t, 0) and ek(t,0) display two crossings, it is 
not difficult to see from the proof of the lemma that no dissipative domains can exist for j 
or k. 

We postpone the proof of this lemma to the end of the section and go on with its conse- 
quences. By applying the results of the previous section we get 

Theorem 6.1 Let H(t, 5) satisfy H4 and H5. If 5 > is small enough, the elements cr(j)j of 
the S -matrix, with a(j) defined in lemma are given in the limit e — > for all j = 1, ■ • • , n, 
by 



ff 1 e~^e~ l k e ^ 5)dz/£ (1 + 0{e)) , a(j) ( > j (6.9) 

where, for a(j) > j (respectively o(j) < j), Cfc, k = j,---,cr(J) =F 1> denotes a negatively 
(resp. positively) oriented loop based at the origin which encircles the set J r (resp. J r 
corresponding to the avoided crossing between ek(t,5) and ek+i(t,6) (resp. e/-—i(z,5)) att r , 
f^ k ek(z,S)dz denotes the integral along ^ of the analytic continuation ofek(0,5) and 9k(Ck) 
is the corresponding factor defined by ( 3.1$ , see figure 



More accurately, with the notations of section we have the improved formula 



JJ e -^) e -'k e ^ S)dz/£ [l + o(e-^)), a(j){ >J . (6.10) 
k=j ' L < J 

The element cr(l)j, I ^ j, are estimated by 

ee Me CTO ,(oo, 5 )-e CT{i) (oo, 5 ))/ e ^ 1m J^S^/e] ^) | > J (6 .H) 

where h is strictly positive (resp. negative) for o~(j) > j (respectively a(j) < j). 
Remarks: 

i) As the eigenvalues are continuous at the degeneracy points, we have that 

lime J Ck kV ' ' =0, VJfe = l,---.p. 6.12 

ii) The remainders 0(e) depend on 8 but it should be possible to get estimates which are 
valid as both e and 5 tend to zero, in the spirit of Q, | MN | and (Rj. 



iii) This result shows that one off-diagonal element per column of the S'-matrix at least can 
be computed asymptotically. However, it is often possible to get more elements by making 
use of symmetries of the S'-matrix. Moreover, if there exist dissipative domains going above 
or below other eigenvalue degeneracies further away in the complex plane, other elements of 
the S'-matrix can be computed. 

iv) Finally, note that all starred quantities in ( 6.1C| ) depend on e. 



Proof: The first thing to determine is whether the loops are above or below the real axis. 
Since the formulas we deduce from the complex WKB analysis are asymptotic formulas, it 
suffices to choose the case which yields exponential decay of s a ^j. It is readily checked in 
the proof of lemma |6.1| below that if o~(j) > j, Dj is above the real axis and if o~{j) < j, 
Dj is below the real axis. Then it remains to explain how to pass from the loop rjj given 



in proposition 4.1 to the set of loops k = j, - ■ ■ ,o~(j) — 1. We briefly deal with the case 
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Figure 5: The loops rjj and (k, k = j, • • • , a(j) — 1. 



a(j) > j, the other case being similar. It follows from lemma 6.1 that we can deform rjj into 
the set of loops £fc, each associated with one avoided crossing, as described in figure §. Thus 
we have 



E 



for the decay rate and, see |3.10| , 



(6.13) 



(6.14) 



for the prefactors. Let i/j be a negatively oriented loop based at t r which encircles J r as 
described in lemma S.l. Consider now the loop associated with this avoided crossing and 



deform it to the path obtained by going from to t r along the real axis, from t r to t r along 
Vj and back from t r to the origin along the real axis. By the point ii) of the lemma we get 



e j (0,S) = e j+1 (0,S) 



along (j, and, accordingly (see (3.14)) 



^•(0,<5)=e-^Vi+i(0,5). 



(6.15) 



(6.16) 



This justifies the first factor in the formula. By repeating the argument at the next avoided 
crossings, keeping in mind that we get ej+i(0, 5) at the end of (j and so on, we get the final 
result. The estimate on s CT (z)j is obtained by direct application of lemma 6.1. □ 



Proof of lemma |6.1| : In the sequel we shall denote " Jj" by a "/ ". We have to consider the 
analyticity properties of ej(z,5) and define domains in which every point z can be reached 
from — oo by means of a path "y(u), u e] — oo,i], j(t) = z such that Im Ajk("y(u), 5) is non 
decreasing in u for certain indices j ^ k when 5 > is fixed. Note that by Schwarz's principle 
if 7(n) is dissipative for {jk}, then ^(u) is dissipative for {kj}. When j(u) = 71 (u) + ^72(1*) 
is differentiable, saying that *y(u) is dissipative for {jk} is equivalent to 

Re (ej(7(u),<5) - e fc (7(ii), S))^ 2 {u) + Im (ej(j(u),S) -e k (~f{u),5))^i(u) > 
VnG]-cxD,t] (6.17) 
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where " ' " denotes the derivative with respect to u. Moreover, if the eigenvalues are analytic 
in a neighborhood of the real axis, we have the relation in that neighborhood 

Im (eAt + is,8) -e k (t + is,6)) = /"" Re (e'At + is', S) - e' k (t + is' ,5))ds' , (6.18) 

Jo 

which is a consequence of the Cauchy-Riemann identity. We proceed as follows. We construct 
dissipative domains above and below the real axis when 5 = and we show that they remain 
dissipative for the perturbed quantities Aj k (z,S) provided 5 is small enough. We introduce 
some quantities to be used in the construction. Let C r G {1, • • • , n} 2 denote the set of distinct 
couples of indices such that the corresponding eigenvalues experience one crossing at t = t r . 
Similarly, N G {1, ■ ■ ■ , n} 2 denotes the set of couples of indices such that the corresponding 
eigenvalues never cross. 

Let I r = [t r — L, t r + L] G R, r = 1, • • • ,p, with L so small that 

min min inf (e'At, 0) - e'At, 0)) = 4c> 0. (6.19) 

r€{l,---,p} {jk}ec r ,j<ktei r J 

This relation defines the constant c and we also define b by 

min min inf leAt, 0) - e k (t, 0)1 > 4b > (6.20) 
re{i,-,p} {jk}eC r ,j<kteR\i r 

, ^ m r „M M*,0)-e fc (t,0)|>46>0. (6.21) 
We further introduce 

I* = { z = t + is\t G I r , \s\ < a}, r = !,-■■, p. (6.22) 

Then we choose a small enough so that the only points of degeneracy of eigenvalues in S a 
are on the real axis and 

min min inf Re (e'Az,0) - e' k (z,0)) > 2c> (6.23) 

r€{l, •■■,p} {jk}£C r ,j<kzeI? J 

min min inf IRe (eAz, 0) - e k (z, 0))| > 26 > (6.24) 

re{i,-, p }{jk}ec r ,j<kzes a \i? Jy 

min inf IRe (eAz,0) - e k (z,0))\ > 2b > 0. (6.25) 

{jk}£N,j<kzeS a 

That this choice is always possible is a consequence of the analyticity of ej(z, 0) close to the 
real axis and of the fact that we can work essentially in a compact, because of hypothesis 
H4. Let a(t) be integrable on R and such that 



a(t) 

— — > max sup 

2 j<ke{i,-,n} | s |< a 



Re (e'At + is, 0) - e' k (t + is,0)) . (6.26) 



It follows from H4 that such functions exist. 
Let r G {1, ■ ■ ■ ,p} and 72 (u) be a solution of 

^ 2 ( u) = _^H ue }-oo,t r -L] 

7 2 (n)=0 u e]t r -L,t r + L[ (6.27) 

[ 72 (n) = +^M ue[t r + L,oo[ 

with 72 (t r ) > 0. Then 72 (ii) > for any u, since 

/ \ / \ ~ \ U a(u')du' lb -i j--, 

7 2 («J = 7 2 (t r Je J tr-L u Gj — 00, t r — L\ 

72(n) = 72(tr) u £}t r -L,t r + L[ (6.28) 

72 («) = y2{t r )JZ+L Au')du'/b u G [t r + L, oo[ 
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and since a(u) is integrable, the limits 



lim 72 (it) = 72 (±00) 

t— »±OQ 



(6.29) 



exist. Moreover, we can always choose 72 (ir) > sufficiently small so that 7 r (u) = u + 
i"f 2 (u) G 5 a , for any real u. Let us verify that this path is dissipative for all {jk} G C r , j < k. 
For u G] — 00, t r — L], we have, using 



Re (e-,(z,0) - e fc (z,0)) < -26 < 0, Vz G S a fl {z|Re z < t r - L} 



(6.30) 



and 



\Im (ej{t + is, 0) - e k (t + is, 0))\ < \s\ sup Re (e^(t + is', 0) - e' fc (t + is', 0)) ,(6.31) 

s'e[o,s] 



(see (gig)), and the definition ( gjg ) 

Re ( ej ( 7 r (u) , 0) - e fc ( 7 r («) , 0)) 7 2 («) + Im (e, ( 7 r (u) , 0) - e* ( 7 r («) , 0) ) 7 i («) = 
-Re ( ei ( 7 »,0) -e fc ( 7 r («),0))^M^ + Im ( ej ( 7 », 0) - e fc ( 7 », 0)) > 
2 72 (n)a(n) - j 2 (u)a(u)/2 > 7^(11)0(11) > 0. (6.32) 

Similarly, when u > t r + L we get, using 

Re (ej(«,0) -e fe (z,0)) > 26 > 0, Vz G 5 Q n {z|Re z > t r + L}, (6.33) 

Re (e,( 7 »,0) - e fc ( 7 r (n), 0)) 72 (n) + Im (ei( 7 », 0) - e fc ( 7 r (u), 0)) 7 i («) = 
Re M 7 »,0) -e fc ( 7 r (n),0)) 72( ^ a(n) + Im (e,( 7 », 0) - e k {f (u), 0)) > 
2 72 (n)a(n) - 72 («)a(«)/2 > 7a(«)a(u) > 0. (6.34) 

Finally, for s G [i r — L, t r + L], we have with ( |6.23D 

/•72(u) 

Im ( ei ( 7 r (n),0) -e fc (7 r (n),0)) = / Re (e^(i' + is, 0) - e^i' + is, 0)) > 

7 2 (n)2c > 7 2 (u)c > 0. (6.35) 

Thus, 7 r (u) is dissipative for all {jk} G C r , j < k. Note that the last estimate shows that it 
is not possible to find a dissipative path for {jk} G C r , j < k below the real axis. 
Consider now {jk} G N, j < k and let 72 (u) be a solution of 



i.e. 



7a(«) = ~ , 72(0) > 0, u G] - 00, +co[, 



l2 (u)= l2 {Q)e'^ a{u ' )du ' /b . 



(6.36) 



(6.37) 



As above, we have 72 (u) > for any u and we can choose 72(0) > small enough so that 
7 + (u) = u + 172(11) G 5 a for any « e R. Since 

Re (ej(z,0) - e k (z,0)) > -26 Vz G 5 Q (6.38) 

we check by a computation analogous to ( 6.32j ) that 7 + (ii) is dissipative for {jk} G N, j < k. 
Similarly, one verifies that if 72 (u) is the solution of 



72(«) = , , 72(0) < 0, u G] - 00, +oo[ 



(6.39) 
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Figure 6: The dissipative domain D and some dissipative paths. 

with |72(0)| small enough, the path 7~(u) = u + £72 (u) below the real axis is in S a for any 
u E R and is dissipative for {jk} E N, j < k as well. 

Finally, the complex conjugates of these paths yield dissipative paths above and below 
the real axis for {jk} E JV, j > k. 

We now define the dissipative domains by means of their borders. Let j + (u) and j~(u), 
u E R, two dissipative paths in S a defined as above with |7^(0)| sufficiently small so that 
7~ is below 7 + . We set 

D = {z = t + is\0 < -7 2 _ (t) < s < 7^(t), t E R}. (6.40) 

Let z E D, and j E {1, • • • , n} be fixed. By assumption H5, the set Xj of indices k such 
that {jfc} € C r for some r E {1, • • • ,p} consists in values k satisfying j < k or it consists 
in values k satisfying j > k. Let us assume the first alternative takes place. Now for any 
k E {1, • • • , n}, there are three cases. 

1) If k E Xj, then there exists a dissipative path Y E -D for {jfc} E C r , j < k constructed 
as above which links —00 to z. It is enough to select the initial condition 72 (t r ) suitably, see 
figure |(| 

2) Similarly, if j < k $ Xj, there exists a dissipative path j + E D for {jk} constructed as 
above which links —00 to z obtained by a suitable choice of 72(0). 

3) Finally, if k > j, we can take as a dissipative path for {jk}, the path 7" E D constructed 
as above which links -00 to 2 with suitable choice of 72(0). Hence D is dissipative for the 
index j, when 5 = 0. If j is such that the set Xj consists in points k with k > j, a similar 
argument with the complex conjugates of the above paths shows that the domain D below 
the real axis is dissipative for j when 5 = 0. 

Let us show that these domains remain dissipative when 5 > is not too large. We start 
by considering the analyticity properties of the perturbed eigenvalues ej(z,5), 5 > 0. Let 
< a' < a be such that 

if n (D UD) = ® , Mr =!,■■-, p. (6.41) 

The analytic eigenvalues ej(z,0), j E {1, • • • ,n}, are isolated in the spectrum of H(z,0) for 
any z E S a , where 

S a = S a \U r=1 ,.., p lf. (6.42) 

For any j = 1, • • • ,n we get from perturbation theory [iKj], that the analytic continuations 
ej(z,5) of ej(ti — L,5) in S a are all distinct in S a , provided 5 is small enough. This is due 
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to the fact that assumption H5 implies the continuity of H(z,5) in 5 uniformly in z G S a , 
as is easily verified. More precisely, for any fixed index j, the eigenvalue ej(t,5) experiences 
avoided crossings at the points {t n , ■■■ ,t rj }. We can assume without loss of generality that 

4 Q 'n/f'=0, Vfc^eK-.rj}. (6.43) 

Hence, for 5 > small enough, the analytic continuation ej(z,5) is isolated in the spectrum 
of H(z, 5), uniformly in z G S a \ U r=ri) ... )7 .. 7° . Since by assumption H5 there is no crossing 
of eigenvalues on the real axis when 5 > 0, there exists a < g < a' , which depends on 5, 
such that ej(z, 5) is isolated in the spectrum of H(z, 5), uniformly in z G S J a , where 

Si = S a \U r=ri ,.. jr . (J r uX) (6.44) 

and 

J r = if n {z\ Im z > g} , r = l,---,p. (6.45) 

Hence the singularities of ej(z,5) are located in {J r=rit ... trj (J r U J r ), which yields the first 
assertion of the lemma. 

Consider a path v r from t r — L to t r + L which goes above J r , where t r is an avoided 
crossing between ej(t, 5) and efc(i, 5), k = j ± 1. By perturbation theory again, ej(t r — L,<5) 
and efc(t r — L, 5) tend to ey(t r — L, 0) and ek'{t r — L,0) as 5 — > 0, for some j', fc' G 1, ■ • • , n, 
whereas ej(t r + L, 5) and efc(i r + L, 5) tend to e}.i(t r + L, 0) and ej'(t r + L, 0) as 5 — > 0, see 
figure ||]. Now, the analytic continuations of ej(t r — L,5) and ek(t r — L,S) along u r , ej(z,5) 
and ek(z,5) tend to the analytic functions ej>(z,0) = ej>(z,0) and e/%'(z,0) = e k '(z, 0) as 
5 — > 0, for all 2 G f r . Thus, we deduce that for 5 small enough 

ej(t r + L,5) = e k (t r + L,S), (6.46) 

since we know that e.j(t T + L,S) = e a (j) (tr + L, 5), for some permutation a. Hence the point 
iii) of the lemma follows. 

Note that the analytic continuations ej(z, 5) are single valued in S a . Indeed, the analytic 
continuation of ej(t r — L,5) along Xv, denoted by ej(z, 5), Vz G T^, is such that 



ej(t r + L,8) = ej(t r + L,S) = ej(t r + L,S) = e k (t r + L,S), (6.47) 

due to Schwarz's principle. We further require 5 to be sufficiently small so that the following 
estimates are satisfied 

min min inf |Re (ej(z, 5) — e/-(z, <5))| > b > (6.48) 

re{l,-,p} W€Or zeSa \I? 

min inf |Re (e,(z,5) - e k (z,S))\ >b>0 (6.49) 
j<k 2fcJ ° 



max sup 

j<fce{i,-,n} Im ,^ 



Re @j(z,8) - e k (z,5)) < a(Re z). (6.50) 



and, in the compacts I" = /°\/° , 

min min inf |Im (ej(z,5) —e k (z,S))\ > 

r£{l, •■■,»} fi fc }6Cr .era 

— min min inf |Im (ej(z,0) — ej.(z,0))\ > |Im z\c (6.51) 

2 rejl, ■■-,»} {i fc }6C r pro, 

max max sup llm (e<j(z,6) — ei~(z,S))\ < 
re{l,-,p}j<ke{l,-,n} ze Ja 

2 max max sup |Im (ej(z,Q) — e k (z,0))\ < |Im z|a(Re z). (6.52) 
re{i,-,p}j<ke{l,-,n} ze J a 
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The simultaneous requirements (6.26) and ( |6.50| ) is made possible by the continuity properties 
of H'(z,S) and the uniformity in 5 of the decay at ±00 of H(z,S) assumed in H4 together 
with the fact that a(t) can be replaced by a multiple of a(t) if necessary, to satisfy both 
estimates. The condition on 5 is given by the first inequalities in ( |6.51| ) and (6.52), whereas 
the second ones are just recalls. 

Then it remains to check that the paths 7 r , 7 + and 7" defined above satisfy the dissi- 
pativity condition ( |6.17 ) for the corresponding indices. This is not difficult, since the above 
estimates are precisely designed to preserve inequalities such as (|6.32 ), ( |6.34| ) and ( |6.35 ). 
However, it should not be forgotten that in the sets I" the eigenvalues may be singular so 
that ( |6.18| ) cannot be used there. So when checking that a path parameterized as above by 
u € R is dissipative, it is necessary to consider separately the case u E R\(U r= i j ... iP I r ), where 
we proceed as above with ( 6.4S| ), ( |6.49 ), ( 6.50 ) and ( 6.1§| ) and the case u € U r= i j ... iP / r , where 
we use use (6.51) and ( 6.52| ) instead of ( 6.18| ) as follows. If u £ I r for r such that t r is a 
crossing point for ej(t,0) and efc(i,0), one takes (|6.51| ) to estimate Im (ef(z, 5) — e# (z, 5)) for 
the corresponding indices f and k', and if t r is not a crossing point for ej(t,0) and ek(t, 0), 
one uses (|6.52D to estimate Im {ef(z, 6) — ek>(z,5)). Consequently, the domains D and D 
defined above keep the same dissipativity properties when 5 > is small enough. 

Let us finally turn to the determination of the associated permutation a. As noticed 
earlier, the eigenvalues ej(z,S) are continuous in 5, uniformly in z G S a . Hence, since the 
eigenvalues 



^(z,0) are analytic in S a , we have 



lim e, (oo, 8) 



, n. 



ej(oo,0) J = 1,2, 
Whereas we have along the real axis (see figure I 



with a defined in the lemma, from which the result follows. 



(6.53) 

(6.54) 

□ 



7 Applications 



Let us consider the time-dependent Schrodinger equation in the adiabatic limit. The relevant 
equation is then Q2.1| ) where H(t) = H*(t) is the time-dependent self-adjoint hamiltonian. 
Thus we can take J = I in proposition |2.1| to get 

H(t) = H*(t) = H*(t). (7.1) 

The norm of an eigenvector being positive, it remains to impose the gap hypothesis in H3 
to fit in the framework and we deduce that the S'-matrix is unitarity, since R = I. In this 
context, the elements of the S'-matrix describe the transitions between the different levels 
between t = —00 and t = +00 in the adiabatic limit. 

We now specify a little more our concern and consider a three-level system, i.e. H(t) = 
H*(t) £ M^(C). We assume that H(t) satisfies the hypotheses of corollary |6.1| with an extra 
parameter 5 which we omit in the notation and displays two avoided crossings at t\ < t2, as 
shown in figure 0. The corresponding permutation a is given by 

<t(1) = 3, a(2) = l, (7(3) =2. (7.2) 



By corollary |6.1| , we can compute asymptotically the elements S31, S12, S23 and Sjj, j = 1, 2, 3. 
Using the unitarity of the S'-matrix, we can get some more information. Introducing 



ImJ e 3 (z)dz 



J = 1,2, (7.3) 
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Figure 7: The pattern of avoided crossings in the adiabatic context. 



where Q is in the upper half plane, with the notation of section ^, it follows that 

S31 = (e-^+W) , s 12 = 0(e- r ^), S23 = 0(e- r ^) (7.4) 

and 

Sjj = l + 0(s), j = 1,2,3. (7.5) 

Expressing the fact that the first and second columns as well as the second and third rows 
are orthogonal, we deduce 



-^g(l + 0(e- 2r2/£ )) (7.6) 
^(l + 0(e- 2r ^)). (7.7) 



•S32 

S22 

Finally, the estimate in corollary ^O] yields 

s 13 = O ^ ee - | ' l|(e2(oo ' <5) - ei(oo ' 5))/e e- r2/E ) = O (^-^2+v 2 +K)/e^ 
where K > 0, since we have that r.,- — > as 5 — > 0. Hence we get 

/ sn S12 ^ e -(r a +r a +J0/ e ) \ 

S= -sE^(l + 0{e- 2r ^)) S22 s 23 (7.9) 

V S3 i -^||(l + o(e- 2r i/ £ )) s 33 ) 

where all Sjk above can be computed asymptotically up to exponentially small relative error, 
using ( |6. 10 ) . 

The smallest asymptotically computable element s\ 3 describes the transition from ei(— 00, 5) 
to e 3 (+oo,5). The result we get for this element is in agreement with the rule of the thumb 
claiming that the transitions take place locally at the avoided crossings and can be consid- 
ered as independent. Accordingly, we can only estimate the smallest element of all, S13, 
which describes the transition from e 3 (— 00, 5) to e\(+oo,5), for which the avoided crossings 
are not encountered in "right order", as discussed in RHPfl . It is possible however to get an 
asymptotic expression for this element in some cases. When the unperturbed levels e2(z,0) 
and e 3 (z,0) possess a degeneracy point in S a and when there exists a dissipative domain for 
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the index 3 of the unperturbed eigenvalues going above this point, one can convince oneself 
that S13 can be computed asymptotically for S small enough, using the techniques presented 
above. 

Our second application is the study of the semi-classical scattering properties of the 
multichannel stationnary Schrodinger equation with energy above the potential barriers. 
The relevant equation is then 

- e 2 ^(t) + V{t)m = E*(t), (7.10) 

where t S R has the meaning of a space variable, <3?(t) E C m is the wave function, e — > 
denotes Planck's constant, V(t) = V*(t) £ M m (C) is the matrix of potentials and the 
spectral parameter E is kept fixed and large enough so that 

U(t) = E — V(t) > 0. (7.11) 

Introducing 

*«>-($«>) 6 ^ (7 ' 12) 

we cast equation ( |7.10| ) into the equivalent form ( |2.1| ) for ifi(t) with generator 

m = {u(t) ) eM2m(c) - (7 - 13) 

It is readily verified that 

H(t) = J _1 iJ*(i)J (7.14) 

with 

J=(° I Q ). (7.15) 

Concerning the spectrum of H(t), it should be remarked that if the real and positive eigen- 
values of U(t), kj(t), j = 1, • ■ ■ ,m associated with the eigenvectors Uj(t) G C m are assumed 
to be distinct, i.e. 

0<kj(t) <kj(t) <■■■ <k 2 m (t), (7.16) 

then the spectrum of the generator H(t) given by ( |7.13| ) consists in 2m real distinct eigen- 
values 

- k m (t) < -km-i(t) <■ ■ ■ < -fci(t) < h(t) < k 2 (t) < ■ ■ ■ < k m (t) (7.17) 
associated with the 2m eigenvectors 



ff(<)x*(() = ±%(t)x*(i). (7.18) 

We check that 

( Xj ± (0),x J ± (0))j = ±2^(0)11^(0)11 ^0, j = l,...,m (7.19) 
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where ||Mj(0)|| is computed in C m , so that proposition 2.1 applies. Before dealing with its 
consequences, we further explicit the structure of 5. Adopting the notation suggested by 
fl7l7l ) and (ffTTsD we write 



H(t) = em^w-EW*) 



and introduce 



n t n t 



(7.20) 
(7.21) 



3=1 



i=i 



c±(t) 



/<f(t)\ 



G C m . 



(7.22) 



V<4(*)/ 

Hence we have the block structure 
c + (— oo) \ / 5+4- 5 4 



5 



-co) 



5. 



5_ 



-oo) 
-oo) 



c+(+oo) 
c~(+oo) 



(7.23) 



where 5 ffT G M ro (C), <r, r G {+,-}. 



Let us turn to the symmetry properties of 5. We get from (7.19) and proposition 2.1 that 
-l 



S+- 
5__ 



I O 
O -I 



S + - 
5__ 



I O 

O -I 



s* ++ 
-si. 



-5! 

5*. 



In terms of the blocks S ar , this is equivalent to 

C C* C C* T 

— 1_ _j I — 0-\ — 1 

5 + +5* + - 5+-5* _ = O 

O O O — i-O i_ — J.. 



.(7.24) 



(7.25) 
(7.26) 
(7.27) 



The block 5++ describes the transmission coefficients associated with a wave traveling from 

the right and 5 |_ describes the associated reflexion coefficients. Similarly, 5 and 5^ are 

related to the transmission and reflexion coefficients associated with a wave incoming from 



the left. It should be noted that in case of equation ( 7.10 ) another convention is often used 
to define an 5-matrix, see [Fl|, for instance. This gives rise to a different 5-matrix with 
similar interpretation. However it is not difficult to establish a one-to-one correspondence 
between the two definitions. If the matrix of potentials V(t) is real symmetric, we have 
further symmetry in the 5-matrix. 



Lemma 7.1 Let 5 given by ( [7. 2$ ) be the S -matrix associated with (7.1C) under condition 



If we further assume V(t) 
^t.2\ ), w e get 



V(t), then, taking ip, (0) G R' 



2m 



5 4 



5 — , 5-| — — 5 — 



, m, in 



(7.28) 



The corresponding results for the 5-matrix defined in [Fl] are derived in fMNfl . The proof 
of this lemma can be found in appendix. We consider now (FfTlGh the case U(t) = U*(t) = 
U{t) G M2(R), which describes a two-channel Schrodinger equation. We assume that the 
four-level generator H{t) displays three avoided crossings at t\ < t2, two of which take place 
at the same point t±, because of the symmetry of the eigenvalues, as in figure ||. By lemma 
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Figure 8: The pattern of avoided crossings in the semiclassical context. 



7.1, it is enough to consider the blocks S ++ and S+— The transitions corresponding to 
elements of these blocks which we can compute asymptotically are from level 1 + to level 2 + 
and from level 2~ to level 1 + . They correspond to elements and sf^~ respectively. With 
the notations 



r. 



Im / ki(z)dz 




3 



1,2, 



where Q is in the upper half plane, we have the estimates 



'21 



.7 



1,2. 



It follows from (7.26) and lemma 7.1 that the matrix S++S+_ is symmetric. Hence 



'11 °21 



+ s 



12 a 22 



'21 °11 



+ S 



22 °12 ' 



whereas we get from ( 7.25 ) 



+-„+- 



S ll s 21 + s 12 s 22 — s ll s 21 + S 12 °22 

The only useful estimate we get with corollary 

s +- = o (e-( ri+r2+/f )/ £ ) , K > 0, 
which yields together with ( |7.30| ) in ( |7.31| ) 

Thus, from ( ffT32| and ( |OoD for sff , 



is 



'12 



'21 



'22 



with 



(7.29) 

(7.30) 
(7.31) 
(7.32) 
(7.33) 
(7.34) 

(7.35) 



< K < min(ri, T2). 



(7.36) 
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Summarizing, we have 



S'4 



'11 



>21 



-4i + ^(l + 0(e- K/£ )) 



and 



5 4 



( 

where all elements s CTT 



(e 



6 22 



s 12 

-(ri+r 2 +K)/ e x 



(7.37) 



(7.38) 



> 7 - fc can be asymptotically computed up to exponentially small relative 
corrections using ( 6.10 ). We get no information on the first column of S-\ but the estimate 



( 5.40| ) where necessarily, ( [7.36 ) holds. However, if there exists one or several other dissipative 
domains for certain indices, it is then possible to get asymptotic formulas for the estimated 
terms. 



A Proof of proposition |2Tl| 

A direct consequence of the property 

H{t) = H*{t) = J- 1 H*{t)J (A.l) 



is the relation a(H(t)) = cr(H(t)). Thus, if a(H(0)) C R, then a(H(t)) C R, for all t € R, 
since the analytic eigenvalues are assumed to be distinct and nondegenerate for all t £ R. 
Let ej(0) be the eigenvalue associated with tpj(0). Then, due to the property H(0) = H#(0) 



{<Pj(0),H{0) Vk {0))j = e fc (0)(^(0),^(0))j = ei (0)(^(0),^(0))j, (A.2) 

for any j,k = 1, • • • ,n. For j = k we get from the assumption (yjj(0), ipj(0))j ^ that 
ej(0) £ R and from the fact that the eigenvalues of H(0) are distinct 

(^•(0),^(0))j = 0, j + k. (A.3) 

The resulting reality of ej(t) for all t £ R and j = 1, • • • , n yields together with ( |A.1| ) 

p^t) = j-'p+m (a.4) 

Hence, using the fact that the P* are projectors, 

n n n 

i=i i=i j'=i 

n 

= -j- 1 ^p;(t)p;'(t)j = -j~ l K*{t)j. (A.5) 

3=1 



Let * £ C n and W(t) be defined by (see (gj)) 

W'(t) = K(t)W(t), W(0) = I. (A.6) 
Then we have 

(W(t)$, W(t)^)'j = (W'(t)$\ JW(t)*) + (W(t)$\JW'(t)tf) 

= (K(t)W(t)<f>\JW(t)V) + (W(t)$\JK(t)W(t)V) 

= {W(t)<f>\J(J- 1 K*(t)J + K(t))W(t)-$) =0. (A.7) 
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Thus, in the indefinite metric, the scalar products of the eigenvectors of H(t), (fj(t) 



W(t)<pj(0) (see (]3.7D ), are constants 

(^•(t), V 9 fc (t)) J = (^(0),^(0)) J . (A. 
We can then normalize the <fj(0) in such a way that 

(<Pj(t), ipk(t))j = (<Pj(0), <Pk(0))j = SjkPj (A. 



with pj £ {+1, —1}- Let i(j(t) and x(t) be two solutions of ( |2.ip . By a argument similar to 
the one above using (|A.1|) , we deduce 

(x(t),m)j = (x(o),m)j- (a.io) 

Inserting the decompositions 

n r* 

m = Y, c ^ iJ ° eAt)dt/£ ^it) (A.H) 

X {t) = Y,dj{t>- i! ° e > {t)M/e ?M (A.12) 

3=1 

in this last identity yields 

3,k=l 3 
n n 

= Y,dj(P)pM0) = E^( ±0 °)^ c i( ±00 )- (A.13) 

3=1 3=1 

Since the initial conditions for the coefficients 

Cj-(-oo) = S jk , d,-(-oo) = (A.14) 

imply 

Cj(+oo) = Sjfc , dj(+oo) = sji, (A.15) 
we get from ( |A.13| ), introducing the matrix R = diag (pi, p2, ■ ■ • , p n ) £ -^n(C), 

= (A.16) 
which is equivalent to the assertion S 1-1 = RS*R. □ 

B Proof of lemma [771] 

Let G = G* = G" 1 be given in block structure by 

G =(l ° T )eM 2m (C) (B.l) 



o -i 

and H(t) be given by with L/"(t) = = 17* (t). Since 



GH{t)G = -H{t) , H(t) = H(t) , (B.2) 
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and the eigenvalues of H(t) are real, it is readily verified that 



GPf{t)G = Pf{t), Pf(t) = Pf{t),j = l,---,m. (B.3) 

Hence 



K{t) = £ PJ'(t)PJ(t) = K(t) = GK(t)G, (B.4) 



T = ± 



from which follows that the solution W(t) of 

W'(t) = K(t)W(t) , W(0) = I (B.5) 
satisfies 



W(t) = W{t) = GW(t)G. (B.6) 

As the matrix of potentials U(0) is real symmetric, its eigenvectors Uj(0) may be chosen real, 
so that we can assume that 

^=L,^JeR 2m (B.7) 



±Jfe,-(0)uj(0) 



Thus it follows from the foregoing that 

ipf{t) = W(t)tf^(0) G R 2m VteR (B.8) 
and satisfies 

G(p f(t) = GW{t)GOpf{<S) = W(t)G<pf(Q) = <pj(t). (B.9) 
Finally, the main consequence of flB.2j ) is that if ip(t) is a solution of 

ieip'(t) = H(t)ip(t), (B.10) 



then (p(t) = Gip(t) is another solution, as easily verified. Thus we can write with (|7T2l| ) , (p~8| 
and (|R9|) that 



m t m t 

3=1 3=1 
m t m , 

= £ c+WtpjWe'fo k ^ dt '/ £ + cJtytp+Me-tf* k ^ dt '/ £ , (B.ll) 

3=1 3=1 

i.e. 



4(t) = c -(t) 



dj(t) = cj(t), Vj = !,-••, m, Vt G R. (B.12) 



Finally, using the definition (7.23) and the above property for t = ±oo, we get for any 
d ± (-oo) G C m 

f (+») \ (s ++ s_ \ f d+(-oo) w^r ^ ( d + (-=c) x 

from which the result follows. □ 
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